Nonclassicality of states and measurements by breaking classical bounds on statistics 
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We derive exceedingly simple practical procedures revealing the quantum nature of states and 
measurements by the violation of classical upper bounds on the statistics of arbitrary measurements. 
Data analysis is minimum and definite conclusions are obtained without evaluation of moments, or 
any other more sophisticated procedures. These nonclassical tests are independent of other typical 
quantum signatures such as sub-Poissonian statistics, quadrature squeezing, or oscillatory statistics. 
This approach can be equally well applied to very diverse situations such as single- and two-mode 
fields, observables with continuous and discrete spectra, finite- and infinite-dimensional systems, 
and ideal and noisy measurements. 



I. INTRODUCTION 

Nonclassicality is a key concept supporting the neces- 
sity of the quantum theory P, 0, B 13, H, H, H H] • A cus- 
tomary signature of nonclassical behavior is the failure of 
the Glauber-Sudarshan P phase-space representation to 
exhibit all the properties of a classical probability den- 
sity. This occurs when P takes negative values, or when 
it fails to be a proper function becoming a generalized 
function or distribution. 

Within standard quantum theory, quantum states play 
two dissimilar but complementary roles: (i) they express 
the state of the system, and (ii) they determine the statis- 
tics of measurements by projection on the system state, 
such as, for example, photon-number and quadrature 
measurements in quantum optics. We may refer to them 
as measured and measuring states respectively. 

In this work we derive exceedingly simple and robust 
practical procedures to reveal the quantum nature of 
measured and measuring states. In this regard, while 
characterization of nonclassical (measured) states has 
been well developed 0, S 1 0, S , much less attention 
has received the characterization of measurements Q. 
One of the purposes of this work is to contribute to fill 
this gap addressing the characterization of nonclassical 
measurements, i. e., when the measuring state is nonclas- 
sical. More specifically, measurements are described by 
positive operator- valued measures (POVMs) A m , such 
that the statistics of the measurement is p m = tr(A m p), 
where p is the measured state. We will say that the mea- 
surement is nonclassical when the P representative of 
some A m takes negative values or is a generalized func- 
tion. In most practical situations, A m define legitimate 
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measuring states p m oc A m so that the measurement is 
nonclassical if and only if some p m is nonclassical. Non- 
classicality of measurements has been recently related 
with the noncontextuality problem in Ref. Q. 

The main contributions of this work are as follows: 

(i) We derive exceedingly simple practical procedures 
that can reveal the quantum nature of states and mea- 
surements. These are upper bounds on measurement 
statistics which are satisfied by all states and measure- 
ments for which the P representative is a non-negative 
function compatible with classical physics. The lack of 
compliance of these statistical bounds is thus a nonclas- 
sical signature. 

(ii) This approach can be applied to arbitrary mea- 
surements, which may involve for example single- or two- 
mode electromagnetic fields, observables with continu- 
ous or discrete spectrum, systems on finite- or infinite- 
dimensional spaces, ideal or noisy measurements, etc. 
(Some of these possibilities are considered in detail be- 
low.) This is in sharp contrast with other nonclassical cri- 
teria that refer exclusively to specific measuring schemes. 

(iii) A key point of this approach is that data analysis 
is reduced to minimum. At difference with other tests of 
nonclassical behavior, in our case definite conclusions can 
be obtained without evaluation of moments, or any other 
more sophisticated data elaborations [E S S 0> S @1- 
This is reflected on the robustness under practical imper- 
fections that may even favour observation of nonclassical 
behavior. 

(iv) These nonclassical tests are independent of other 
typical quantum signatures of nonclassical behavior such 
as sub-Poissonian statistics, squeezing, or oscillatory 
statistics [T| ■ To this end we propose examples of quan- 
tum states violating classical bounds that present no such 
typical quantum signatures. 

To derive the nonclassical tests we will use the P and 
Q phase-space representatives associated to any operator 
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A, defined as 

A = J d 2 aP(a)\a){a\, Q{a) = ^(a\A\a) , (1.1) 

where \a) are coherent states, a\a) = a\a), and a is the 
annihilation or complex-amplitude operator. They are 
suitably normalized 

J d 2 aP(a) = J d 2 aQ(a) = trA, (1.2) 

with d 2 a = dxdy, where x, and y are the real and imag- 
inary parts of a = x + iy. The measured statistics 
p m = tr(A m p) can be then expressed as 

Pm = 7T J d 2 aP m (a)Q{a) = it J d 2 aP(a)Q m (a), 

(1-3) 

where P{a) and Q{a) are the P and Q representatives 
of the measured state p, while P m (a) and (J m (a) are the 
ones associated to the POVM A m . 

In Sees. II and III we derive simple bounds to p m 
able to reveal the nonclassical nature of measuring states 
p m oc A m and measured states p, respectively. The ro- 
bustness of these criteria under practical imperfections 
is examined in Sec. IV. This formalism is further ex- 
tended to two-mode situations in Sec. V, and adapted to 
finite-dimensional systems in Sec. VI. 



II. NONCLASSICAL MEASUREMENTS 

From Eq. (|1.3p we can derive classical bounds dis- 
closing nonclassical measurements. For every ordinary 
non-negative function P m (a) > it holds that for every 
a 

P m {a)Q{a) < P m (a)Q max , (2.1) 

where Q max is the maximum of Q(a) (note that Q{ot) is 
always a positive and well behaved function). Applying 
this to the first equality in Eq. (|1 .3[) we get the following 
upper bound for p m , provided that trA m is finite, 

Pm < 7T<2maxtrA m . (2.2) 

Equation (|2.2|) can be violated if P m (a) fails to be posi- 
tive or when it becomes a generalized function. In both 
cases Eq. (|2.1[) fails to be true. Therefore, the violation 
of condition (|2.2|) is a signature of nonclassical measure- 
ment. 

The existence of P m (a) as a classical probability den- 
sity for all m allows us to understand the measurement as 
a classical stochastic process @ between the phase space 
and the sample space, with transition probability kernel 
given by K(m,a) — nP m (a). Conversely, the failure of 
K(m, a) to be a classical conditional probability density 
denotes the quantum nature of the measurement process. 



In order to detect the violation of the classical bound 
(|2.2p the only prior information required about the mea- 
surement being performed is the trace tr(A TO ). This can 
be measured using explicit practical methods (see some 
proposals in the Appendix). In any case this is not a very 
stringent condition since in most practical situations this 
can be inferred from simple rough analyses of the exper- 
imental arrangement, by symmetry considerations, etc. 

Note that coherent states \a) are useless as mea- 
sured states to reveal nonclassical measurements since 
TrQmax = 1 so that Eq. ()2.2() leads to the trivial bound 
Pm < trA m for all measurements (10|. This bound is 
trivial because, using the Cauchy-Schwarz inequality, 

|tr(ylS t )| 2 < ti(AA^)tr(BB^), (2.3) 

we get 

p 2 m = [tr ( P A m )} 2 < tr (p 2 ) tr (A 2 ,) , (2.4) 

and using that for positive operators tr(A 2 ) < (trA) 2 we 
get 

Pm = tr(pA m ) < tip trA m = trA m . (2.5) 

Otherwise, quantum or classical state other than coher- 
ent may be used since the weight of the criteria relies on 
the behavior of P m {a). 

This approach is next illustrated with the examples of 
photon-number and field quadrature measurements per- 
formed on a single-mode electromagnetic field. 

A. Photon-number measurements 

In order to illustrate this formalism the simplest ex- 
ample is the ideal photon- number measurement, A„ = 
p„ = \n)(n\, where \n) are number states, a)a\n) = n\n) 
so that trA„ = 1. In this case the classical bound in Eq. 
(|2.2|) becomes 

Pn < TrQmax = Pb, (2.6) 

which is actually independent of the outcome n. 

A readily demonstration of the nonclassical nature of 
the photon-number measurement is provided when n = 1 
and the measured state is the one-photon state \n = 1). 
In such a case p± = 1, 

Q(a) = 1^- exp(-|a| 2 ), 7rQ max = -, (2.7) 
7r e 

where the maximum occurs for |a| = 1. Thus we have 
that 

Pi = 1 > 7rQ max trAi = -, (2.8) 

e 

so that the measurement is nonclassical and the classical 
upper bound is surpassed by 172 %, since (pi —pb)/pb = 
1.72. 
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B. Quadrature measurements 

Concerning quadrature measurements (implemented 
in practice by homodyne detection [l|) we have A x = 
\x)(x\, where \x) are the eigenstates of the quadrature 
operator 

X = ^(a j +a), X\x)=x\x), (2.9) 

being the optical analog of mechanical position or linear 
momentum. In this case trA x is not finite since \x) are 
not normalizable (x\x') = 8(x — x'). 

In order to avoid this difficulty we can appreciate that 
the P representative of |x)(a;|, P x (a = x' + iy'), docs 
not depend on y' . This is an observable property, for 
example, via the independence of statistics under dis- 
placements of the measured state along this coordinate. 
Thus we can rearrange Eq. (| 1 - 3[) in the form 



Px = tt / dx'P x (x')Q(x'), Q(x') = I dy'Q(x',y'), 

(2.10) 



so that Eq. (|2.2p is replaced by 

Px — ^QniaxX^x^x-j 



(2.11) 



where Q max is the maximum of Q{x) when x is varied, 
and 

tr.A, = J dx'P x (x') = J dx'Q x (x'), (2.12) 
where P x , and Q x are the representatives of A x , with 
Q x (x') = -\(x\a = x' + iy')\ 2 

7T 

= -J -exp \-2(x-x') 2 ] . (2.13) 

7T V 7T 

This leads to tr x A x — l/n and to the classical upper 
bound, 

Px < Qmax = Pb, (2-14) 

that does not depend on the outcome x. 

1. Thermal- chaotic state 



In order to look for violations of bound (|2.14|1 let us 
consider that the measured state p is the thermal-chaotic 
state whose expression in photon-number basis is 



p tc = (i-o£e»R 



(2.15) 



n=0 



where £ is a real parameter with < £ < 1. These states 
describe most classical light sources. The mean number 
of photons n tc and the quadrature variance are 

ntc = j^, (AX) 2 = 1 (1 + 2n tc ) , (2.16) 



while the Q and Q functions are 

= ¥(^+T) ex P (-£Ei) 



Q{x') = 



■■ exp 



so that the upper bound in Eq. (|2. 14(1 reads as 

1 

\Jix[nt c + 1) 



(2.17) 



(2.18) 



The statistics of the quadrature measurement p x 
|(a;|p tc |a;)| 2 is Gaussian 



1 



Px 



2irAX 



exp 



2{AX) 2 



(2.19) 



and the output most likely to break bound (12. 2| is x — 0, 
since it maximizes p x . This outcome will infringe the 
bound provided that 



1 



Po 



■ Pb 



1 



yjn{ntc + \) V*(. n tc + 1) 



(2.20) 



which holds for every ntc- In particular for nt c — (the 
vacuum state) we have po = 0.80 and pb — 0.56, so that 
the classical upper bound is very clearly surpassed by 
I00( j5o -Ph)/P6 = 43%. 

The outputs x that contravene Eq. (|2.14|) are all x 
such that 



x 2 < (AX) 2 In 



1 



4 (AX) 



(2.21) 



For n tc = these are all x in the interval —0.42 < 
x < 0.42, which occur with a 60 % probability since 

I-0A2P xdx - °- 60 - 

2. Squeezed vacuum 

As a further example, when the measured state is the 
squeezed vacuum the quadrature statistics has again the 
Gaussian form (|2.19[) . being the Q function 



Q(x,y) = --T— -rTX-^^i' 



7T 1 +4(AX) 2 



2x 2 + 8(AX) 2 y 2 



so that 

Q(x) = 
and 



tt[1 + A{Axy 



■ exp 



1 + 4(AI) 5 



2x 2 
1 + 4(A1) 5 



tt[1 + 4(AX) 2 ]' 



(2.22) 
(2.23) 

(2.24) 
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The output most likely to break bound (|2. 14[) is x = 0, 
and in such a case the classical bound is surpassed for all 
AX, since 

p ° = T^Kx > Pb = ^[i + 4(ax)T (2 ' 25) 

The most favorable situation is when AX is as small as 
possible. For example, for AX = 0.1 we have po = 4.0 
and Qmax = 0.8, so that there is a percentage of violation 
of 100(po— Pb)/Pb — 400% approximately. The outputs x 
that contravene Eq. (|2.14p are given by Eq. (|2.2ip , which 
for AX = 0.1 is the interval -0.18 < x < 0.18, that 
represents the 93 % of all outcomes since f_ 01 ~p x dx ~ 
0.93. 




2 4 6 8 10 

n 



FIG. 1: Classical upper bound (|3.5|l for the probability of 
detecting n photons on a classical state. 



A. Photon-number measurement 



III. NONCLASSICAL STATES 

In this section we derive classical bounds disclosing 
nonclassical measured states. They can be derived from 
the last equality in Eq. (|1.3p by considering that for 
classical states, i. e., for ordinary non- negative functions 
P(a) > 0, we get 

P(a)Q m (a) < P{a)Q 

m , max j 

(3.1) 

where Q m ,max is the maximum of Q m (a). Applying this 
to the last equality in Eq. (|1.3[) and taking into account 
Eq. (|1.2[) , we get the following upper bound for p m , 

Pvri — ^Qm.raaxj (3-2) 

that holds for every P(a) compatible with classical 
physics. If this condition is violated for any m the state 
is not classical. The bound becomes an equality when 
the measured state is the coherent state |a m . max ) with 

Qm.max — Qm(^m,max)- 

Moreover, when the measured state is pure p = \ip) (ip\ 
we have that \ip) is nonclassical if and only if there is at 
least a measurement for which the classical bound (|3.2p 
is violated. The violation of bound (|3.2p is clearly a suffi- 
cient condition. This is also necessary since for every non- 
classical \ip) we can consider a POVM with A = 
In such a case the upper bound (|3.2[) is surpassed because 
Po = 1 while 7rQo,max = | (a max | V>) | 2 < 1, since otherwise 
the equality |(a max |^>)| = 1 would imply that l^) is a 
coherent state and thus classical. 

Note that the POVM A Q = \a){a\/ir defined by the 
coherent states \a) (implemented in practice by double 
homodyne and heterodyne detection) are useless for the 
detection of nonclassical states, since TrQ a ,m&x = 1 and 
bound (|3.2p becomes trivial p a < 1. 

For the sake of illustration we particularize this ap- 
proach to two meaningful practical situations. These are 
photon-number (Sec. Ill A) and quadrature measure- 
ments (Sec. Ill B). Then we apply them to different 
measured states (Sec. Ill C). 



For photon- number measurements A„ = the Q 

function is 

1 \a\ 2n 
Q„(a) = -exp(-H 2 )^-, (3.3) 
7r n\ 

and the maximum occurs at \a\ = \pri: 

1 n" 
Qn.max = - exp(-n) — . (3.4) 
7r n\ 

If the measured state is classical, the photon-number 
statistics p n is thus bounded by 

n n 

Pn < exp(-n)— = p b n . (3.5) 
n\ 

This bound was previously derived in Ref . ■ The upper 
bound p^ n is the probability of detecting n photons in the 
coherent state \a) with \a\ — ^/n, which is the classical 
state for which p n is maximum. In Fig. 1 we have rep- 
resented pb, n as a function of n showing that for large n 
it decays as pb, n — l/V2wn approximately, in agreement 
with the Stirling approximation n\ ~ \phrnn n exp(— n). 

1. Independence of sub-Poissonian statistics 

We can show that the nonclassical criterion on photon- 
number measurements Q3.5P is independent of sub- 
Poissonian statistics. The deviation from Poissonian 
statistics is usually assessed by the Mandel parameter 

Q 

Qm = — n 1- 3.6) 

(n) 

The independence holds because: (i) there are sub- 
Poissonian states that satisfy the classical bounds (|3.5|) 
for all n, and (ii) there are super-Poissonian states that 
infringe them. To show this let us consider the state in 
the number basis, 

p=(l-p)\Q)(Q\+p\N)(N\, (3.7) 
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with < p < 1 so that 

(ri)=pN, (An) 2 = N 2 p(l -p), (3.8) 

and 



! M 



N(l-p) - 1. 



(3.9) 



(i) For N = 1 this state is sub-Poissonian for all p with 
Qm = — P < 0, and satisfies the classical upper bounds 
([37g| for all n when p < 1/e. (ii) For p < (N-l)/N the 
state is super-Poissonian since Qm > 0, and infringes 
bound (|3.5p when p > pb.^r. These two requirements are 
compatible since it holds that (N — l)/N > p^N for all 
N > 1. 



B. Quadrature measurement 

As a further example we may consider the measure- 
ment of the quadrature X in Eq. (|2.9p . so that A x = 
\x)(x\. In such a case, from Eq. (|2.13p we get nQ xmax — 
\J2/~k, and the classical upper bound for the statistics p x 
of the quadrature measurement is 



Px<\/- 



Pb, 



(3.10) 



that does not depend on the output x. The maximum 
is obtained for a coherent state \a) with (a + a*)/2 = 
x. For states with Gaussian p x the infringement of Eq. 
(|3.10P is equivalent to squeezing of quadrature X since 
the maximum of p x is l/(AX^/2n) and 



1 



AXV2n 



> 



— AX < i = AX vacuum . (3.11) 

TT Z 



For non Gaussian the situation can be different as 
shown below. 

This example is interesting since quadrature mea- 
surements are more experimentally feasible than num- 
ber measurements. For a further discussion about the 
quantum-classical relation in terms of quadrature distri- 
butions see Ref. fill. 



where < p < 1, leading to a photon- number statistics 



Pn 



•p(i-or + (i-p)*f. 



(3.13) 



For example, for p = 0.5, n = 1, and n tc — 9 the proba- 
bility of detecting a single photon is p\ = 0.545, while the 
upper bound in Eq. (|3.5p for n = 1 is pb.i = 1/e = 0.368, 
so we have a clear infringement of the classical condition 
by 100(pi - p b ,i)/pb,i = 48%. 
In this case the photon-number distribution (|3.13p is 
highly super-Poissonian with Qm — 11.2. Furthermore, 
we can easily show that there is no quadrature squeezing 
since for quadrature operators 



X e = - [a^ exp(—i8) + aexp(i6*)] , 
we have in state (|3.12p that (Xg) — and 



(3.14) 



(AX e 



[pn tc + (l-p)n Q } + -. (3.15) 



For the above parameters, no = 1, p = 0.5, and ntc — 9, 
we get (AXg) 2 = 11/4 for all 6, which is far above the 
upper limit for squeezing (AXg) 2 ^ cuum = 1/4. Finally, it 
can be appreciated that there are no oscillations in the 
photon-number distribution. 

For state (|3.12p the origin of nonclassical behavior is 
that P(a) is always more singular than a delta function 
for all p 7^ 1. This is because 



where 



and 



P(a)=pP tc (a) + (l-p)P no (a), 



p , , i /Id 

Ptc(a) = exp 

TT^tc V n tc 



(3.16) 



hr- > ( 3 - 17 ) 



^(a)={l + l^)s^(a). (3.18) 

Thus P(a) is more singular than a delta function since 
otherwise we would be able to express P na =i( a ) as a li n " 
ear combination of two ordinary functions. 



C. Examples 

Let us consider some meaningful simple examples of 
states violating the classical upper bounds (|3.5p and 

pro]) . 

1. Incoherent superposition of thermal and number 

As a feasible state that can infringe Eq. (|3.5p ^ 
us consider the incoherent superposition of the thermal- 
chaotic state in Eq. (|2.15p and the photon number state 

K) 

p = P(Hc + (l-p)\no){no\, (3.12) 



2. Photon-added thermal state 

The previous states are associated with P representa- 
tives more singular than a delta function. Next we con- 
sider states with nonsingular P{a) function taking neg- 
ative values. This is the case of the single-photon-added 
thermal states that, in the photon-number basis, read as 

1 

OO 

pi = (1 - dWptca = (1 - £) 2 J2 f-^WR (3.19) 

n=l 

where pt c , £ and nt c are in Eqs. (|2.15p and (|2.16p respec- 
tively. 
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The P representative is well behaved but nonpositive, 
1 



violation classical bound 



P(a) 



[(n tc + 1) \a\ 2 - n tc ] exp (-\a\ 2 /n tc ) 



and the photon-number statistics is 

Pn - (1 - 2 C n_1 » 



(3.20) 



(3.21) 



For example, for n = 1,2 the classical upper bounds 
are surpassed provided that 

Pi > Pb,i — 1/e <-> n tc < - 1 = 0.65, 
P2 > Pb,2 = 2/e 2 ^> 0.30 < n tc < 0.82. (3.22) 

Let us show that this nonclassical behavior is indepen- 
dent of other nonclassical features. There is no quadra- 
ture squeezing since for the rotated quadrature operators 
X e in Eq. (|3~14l) we have (AX e ) 2 = (3 + 4n t c )/4 > 1/4 
for all 6. Also, it can be appreciated in Eq. (|3.2ip that 
there is no photon- number oscillations. Finally, the Man- 
del parameter is 



Qm = 



He 



1 



2n tc 



1 



(3.23) 



so we get super- Poissonian statistics for all n tc > 1 /y/2 = 
0.71. 

Therefore, the states with 0.71 < n tc < 0.82 are 
nonclassical since p2 > Pb,2, although they have super- 
Poissonian statistics, present no squeezing, and have no 
oscillatory statistics. 



3. Coherent superposition of coherent states 

Another interesting example is provided by the coher- 
ent superposition of two coherent states with opposed 
complex amplitude [12| (referred to as even and odd su- 
perpositions [Hj]) 



\a±)=N±{\a)±\-a)), 



with 



cx P (H 2 /2) 
2^cosh(H 2 ) ' 



N- 



exp(|a| 2 /2) 
2^sinh(|a| 2 )' 



(3.24) 



(3.25) 



In this case the P{a) is a distribution involving an in- 
finite number of derivatives of the delta function, since 
the normally-ordered characteristic function is a real ex- 
ponential. For definiteness let us focus just on the even 
states \a + ). 

For the even case \ot+) we have the following photon- 
number statistics 



Pn 



, ^ a l/, r >s for n even, 
for n odd. 



(3.26) 




FIG. 2: Percentage 100(po ~ Pb)/Pb of maximum violation 
of the classical bound (|3.10[1 (dashed line) and percentage of 
squeezing 100(1 - 2AX 9 , min ) in Eq. (pQ2)l (solid line) for the 
even state (|3.24[1 as functions of \a\. 



Numerically it can be easily seen that this is inconsistent 
with the upper bound in Eq. (|3 . 5|) for \a\ > 0.64. More- 
over, for |cv| >> 1 we have cosh(|a| 2 ) ~ exp(|a!| 2 )/2 and 
(always for even n) 



■exp(-|a| 2 ), 



(3.27) 



so that for \a\ 2 — n we get that p n is twice the upper 



bound pb /n in Eq. (|3.5[) . 

Concerning quadrature measurements, let us consider 
states with purely imaginary complex amplitude a = 
±i|a| that have the following quadrature statistics: 



Px 



4A 2 



— cos 2 (2|a|a;) exp i 



-2x 2 



(3.28) 



For every \a\ the maximum of p x holds for x — being 
Po = 4N+y/2/n. In Fig. 2 (dashed line) we have rep- 
resented the relative amount of violation of Eq. (|3.10|) . 
100(p — Pb)/Pb, as a function of \a\, showing a 100 % 
violation for large \a\. This is because for \a\ >> 1 we 
have Nl ~ 1/2 so that pa ~ 2^/2/7r, which is twice the 
classical upper bound pb — \J~2pn. In Fig. 2 this is also 
compared with the percentage of squeezing in the same 
state (solid line). 

Next we show that the even states infringe classical 
bounds with super-Poissonian photon-number statistics 
and with negligible quadrature squeezing. Concerning 
photon-number statistics we have 



! tanh(|a| 2 ), 
|a| 2 tanh(|a| 2 ), 



so that 



2M : 



sinh(2|a| 2 ) 



(3.29) 



(3.30) 



and these states are always super-Poissonian (unless a 
0). 
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Concerning quadrature squeezing, the minimum un- 
certainty for rotated quadratures (|3.14[) in the state \a + ) 
when 9 is varied is 



(AX t 



?■ 

/mm 



- (2|a| 2 tanh|a| 2 



2\a\ 



1 



(3.31) 



The percentage of squeezing defined as 



100- 



AX 



vacuum 



AX, 



AX 



100 (1 - 2AX e 



vacuum 



(3.32) 

is represented in Fig. 2 as a function of |a|. These states 
present squeezing only for small |a| being negligible for 
|a| > 2. For instance, for \a\ = 3 we have (AXg tm i n ) 2 = 
0.24999986, which means a fully negligible 2.75 x 10~ 5 % 
squeezing, while Eqs. (|3.5p and (|3.10p are infringed by a 
100% for the same state. 

At difference with the preceding examples in this case 
the 100 % violation of classical bounds for large |a| has 
a simple explanation in terms of the oscillatory char- 
acter of the statistics (|3.26[) and (|3.28p . For large \a\ 
the number and quadrature statistics are the same of co- 
herent states (that would saturate the classical bounds) 
but maximally modulated. Because of normalization, 
the vanishing terms must be compensated by nonvanish- 
ing terms reaching twice the coherent-state values. This 
factor of 2 leads to the 100% violation of the classical 
bounds. 

Let us note that the criteria presented in this work re- 
veal the nonclassical nature of these states for all a, but 
specially clearly for large |a|. This is sharp contrast with 
sub-Poissonian number statistics and quadrature squeez- 
ing, that hold only for small \a\, as illustrated in Fig. 2 
for example. 



IV. EFFECT OF IMPERFECTIONS 

One of the key features of this approach is that the data 
analysis is reduced to minimum. This favours obtaining 
reliable results from non ideal measurements affected by 
imperfections, such as damping, finite efficiencies, or fi- 
nite sampling. We stress that this approach applies to 
any measurement, both ideal and imperfect, so that ex- 
perimental imperfections can be always embodied into 
the measuring POVM. Nevertheless, since imperfections 
usually deteriorate nonclassical properties it is reason- 
able to investigate their effect on the above nonclassical 
criteria. 



A. Inefficient detection 

For definiteness we consider real detectors affected by 
field damping (with bath at zero temperature) and finite 
quantum efficiency, which can be modeled by placing a 
beam splitter of amplitude-transmission coefficient t = 
yfrj in front of a perfect detector, where r\ < 1 represents 
both losses and efficiencies (HI XYm . 




0.2 0.4 0.6 0.8 1 



ri 



FIG. 3: pt,i in Eq. (|4.3|l as a function of r\ for n tc — 0.7. 



1. Nonclassical states 

The effect of the beam splitter for the detection of non- 
classical states can be easily accounted for by computing 
the measured state after the beam splitter pt as 



Pi 



d 2 aP{a)\ta)(ta\ 



(4.1) 



where P(a) is the P function of the measured state. The 
measured statistics becomes 



pt, m = tr (A m p t ) = tt J <T aP(a)Q m (ta) . (4.2) 

Since the maximum of Q m (ta) when a is varied is the 
same as the maximum of Q m (a), there is no change in 
the classical upper bound in the right-hand side of Eq. 
(Ell). 

Nevertheless, imperfections affect the statistics replac- 
ing p m by pt.m m the left-hand side of Eq. (|3.2p . This 
can be easily seen, for example, for photon-added ther- 
mal states p,19p . For inefficient detection the photon- 
number statistics of ideal case (|3.2ip for n = 1 is replaced 
by 



PtA = V- 



1 + 2n t 



V n tc 



(1 + rjn tc y 



(4.3) 



In Fig. 3 we have represented ptA as a function of r\ 
for rite = 0.7. Decreasing r\ from r\ = 1 increases pt,i, 
leading to break the classical bound (]3.5[) in the interval 
0.30 < rj < 0.89. 

It is worth pointing out that, at difference with other 
nonclassical tests, where imperfections degrade nonclassi- 
cal behavior [l5| , in this case larger losses and decreasing 
efficiencies may favour the observation of nonclassical be- 
havior of measured states. This noticeable effect arises 
because imperfections rearrange the probability distribu- 
tion p n , so that with increasing imperfection some prob- 
abilities may increase beyond the classical bounds, as is 
the case of ptA in this example. 
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2. Nonclassical measurements 



so that 



For nonclassical measurements we can follow two dif- 
ferent strategies: (i) we can address the nonclassical be- 
havior of the ideal POVM A m associated to 77 = 1. This 
is the analog of the preceding subsection where we in- 
vestigated the nonclassical properties of the input state 
before being affected by imperfections, (ii) Alternatively, 
as mentioned above we can examine the nonclassical be- 
havior of the effective POVM A m embodying all imper- 
fections as part of the measuring scheme. 

(i) Concerning the nonclassical behavior of the ideal 
POVM A m we can compute the effect of inefficiencies as 

Pt,» = - / cPacP/3 s (a\ a (P\Tp® p T^\l3) a \a) s P m {a), 

(4.4) 

where \a) s are coherent states in the signal mode (the 
mode of the measured state p) and \(3) a are coherent 
states in the auxiliary mode (the other input port of the 
beam splitter assumed in the vacuum state po), P m (a) is 
the P representative of the ideal measurement A m , and 
T is the unitary transformation describing the effect of 
the beam splitter, with 



T^\P) a \a) s = \t(3 + ra) a \ta - rf3) 6 



(4.5) 



being r = \/l — t 2 . This leads to the following form for 
the statistics 



Pi,, 



d 2 aQ(a)P m (a), 



(4.6) 



where Q(a) is defined here as 

Q{a) = [ d 2 PQ (t(3 + ra)Q(ta - r/3), (4.7) 



with Qo and Q being the Q representatives of po and p, 
respectively. Note that in the ideal case 77 = 1 (t = 1, 
r = 0) Q(a) is the Q function of the measured state 
Q(a) = Q(a). From Eq. (|4.6[) we can derive the classical 
upper bound 



Pt,m — 7T^t?maxtrA m , 



(4.8) 



which holds for classical measurements with P m {<^) > 0. 
We can appreciate that finite quantum efficiencies modify 
the classical upper bounds in comparison with the ideal 
detection in Eq. (|2.2[) by replacing Q max by Q ma x- 

Let us illustrate this analysis with the example where 
the ideal POVM is one-photon detection, Ai = |1)(1| 
in the number basis, and the measured state is the one- 
photon state The case 77 = 1 was considered in Sec. 
II A above. When 77 < 1 we have 



p t ,1 = 77, trAi = 1, 



(4.9) 



and 



Q(a) = -[ V (M 2 - 1) + 1] exp (-H 2 ) , (4.10) 



exp 



2?7- 1 
77 



(4.11) 



and the violation of the classical upper bound (|4.8|) oc- 
curs provided that 



exp 



277 — 1 



< 1 



(4.12) 



Therefore, in this example the nonclassical behavior of 
the ideal measurement is disclosed provided that the 
quantum efficiencies are above 50 %. 

(ii) Alternatively, if we embody decaying mechanisms 
and inefficiencies in the effective POVM A m we get from 
Eq. 



Pt,m ^ ^QmaxtrA m , 



(4.13) 



where Q m ax is the maximum of the Q function of the 
measured state p. We can compute trA m taking into 
account the effect of the beam splitter as in Eq. (|4.ip : 

trA m = i J d 2 a(a\A m \a) 

= - / d 2 a(ta\A m \ta) = ^trA m , (4.14) 



leading to 



r. 



n 



(4.15) 



We can appreciate that the effect of imperfections is sim- 
ply expressed by increasing the classical upper bound by 
a factor of I/77. 

Let us illustrate this approach with the same example 
of inefficient one-photon detection, (ideal POVM Ai = 
|1)(1| and a one-photon state |1)) so that after Eq. (|2.7[) 
the classical bound (|4.15p is 



Pt,i = V 



1 

= n < — . 



(4.16) 



Thus the effective POVM A m shows nonclassical behav- 
ior when 7/ > l/\/e = 0.6065. 

We can appreciate that the two approaches (i) and 
(ii) lead to two different bounds, (|4.8[) and (|4. 1 5[> . as 
clearly illustrated by the example of one-photon detec- 
tion. We stress that this difference is natural since the 
classical bound (|4.8|) is sensitive to the nonclassicality 
of the P representative of the ideal POVM A m , while 
bound (|4.15p is sensitive to the nonclassical character of 
the P representative of the effective POVM A m . 



B. Finite sampling 

When the number of measurements N is finite, the 
probability p m becomes an statistical variable that can 
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be expressed as p m {N) = k/N, where the integer k is 
the number of outcomes m after N trials. (This anal- 
ysis applies both to detection of nonclassical states and 
measurements.) The dichotomic character of the mea- 
surement (outcome m with probability p m and outcome 
not to with probability 1 — p m ) implies that k follows the 
binomial distribution 



V k (N) 



N 



p k m {l- Pm ) 



N-k 



(4.17) 



so we have 



(p m (N)) = \f =p. 



A Pm (N) 



(fc) 
JV 

Afc _ j p m {l ~Pm) 

N ~ V N 



(4.18) 



For all the above examples we have roughly p m ~ 0.5, so 
that for N ~ 100 we have Ap m /p m ~ 0.1. Thus, even 
for moderate number of trials, the uncertainty caused by 
finite sampling is clearly below the amount of violation 
of classical upper bounds Pb, m , since (p m — Pb,m)/Pb,m 
is at least five time larger than Ap m /p rn in the above 
examples. 



V. TWO-MODE OBSERVABLES 

The above single-mode approach in Eq. (| 1 . 3[) can be 
easily generalized to two-mode observables by expressing 
the statistics p m = tr(pA m ) as 

2 / j2 .2 



and 



p m = n 2 J d 2 ad 2 pP rn (a, [3)Q{a, (3), (5.2) 

where P, and Q are the two-mode phase-space represen- 
tatives for the measured state p 



Pm = w z / d 2 ad z /3P(a, (3)Q m {a, /3), (5.1) 



p= d 2 ad 2 f3P(a,/3)\a,(3)(a,(3\ 



Q(a,(3) 



t,(3\p\a,/3), 



(5.3) 
(5.4) 



P m , and Q rn refer to the corresponding representatives 
of the POVM A m , and \a,/3) are two-mode coherent 
states. The parameter m represents all the indices neces- 
sary to label the outcomes. The maxima of Q(a,j3) and 
Q m {oL, (3) provide suitable upper bounds for the statistics 
of classical measurements and states, respectively. 

In this regard we note that for bipartite systems non- 
classical P(a,(3) is a necessary condition for entangle- 
ment [16]. For definiteness we focus on the nonclassical 
behavior of states. The analysis of nonclassical measure- 
ments would be analogous. 



A. Nonclassical states by photon-number detection 

For the case of joint two-mode photon-number detec- 
tion we get that for classical states the joint probability 
of detecting ni and 7i 2 photons is simply bounded by the 
product of the one-mode upper bounds 



Pm,n 2 < Pb.mPb.n^ = CXp[-(77i + 77 2 )] 



77l!?72! 



(5.5) 



The maximum for fixed 77! + 77 2 occurs when n± — or 
772 = while the minimum occurs for coincident outputs 
m = 77 2 . 

For the total number 77 = 771 + 772 the statistics is given 
by 

/n 
d 2 ad 2 l3P(a, 13) £ Q m . n - m (a, 0), (5.6) 



m=0 



where 



?nx,n a (a,/9) = 4|(a|77 1 )| 2 K/3|77 2 )| 2 , (5.7) 



^1,2) being number states in the corresponding modes. 
It can be easily seen that 

71 1 I |2n 

V Q m ,n-m(a,0) = — ^Vexp(-| 7 | 2 ), (5.8) 

z — ' 7T Z 77! 

m=0 

with I7I 2 = \a\ 2 + \(3\ 2 . For fixed 77 the maximum occurs 
for I7I 2 = 77 so that 



El 77" 
Qm,n-m{0L, (3) < ~r exp(-n) 
7r n,\ 



m=0 



7T 77! 



and 



Pn < — rexp(-77), 



(5.9) 



(5.10) 



which is equal to the single-mode counterpart (|3.5 



B. Nonclassical states by quadrature-difference 
measurement 

Let us consider the measurement of the quadrature 
difference X = X\ — X2, where -Xj. 2 represent the same 
quadrature operator in each mode, which is described by 
the POVM 



A x = J dx'\x + x')i{x + x'\® \x') 2 (x'\, (5.11) 

where \x)j are the eigenstates of Xj, with j = 1,2. In 
this case we have 



1 



TT 2 y/TT 



exp [-(x- xi +x 2 ) 2 ] , (5.12) 
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where Xj is the real part of ay , so that 



VI. SPIN SYSTEMS 



and the classical bound is 



1 



Px < — = = Pb- 

V7T 



(5.13) 



(5.14) 



Note that the two-mode bound does not depend on the 
outcome x, being lower than the single-mode counterpart 
(|3.10p . For states with Gaussian p x the violation of this 
bound is equivalent to AX < l/y/2. This is equivalent 
to two-mode squeezing since for pairs of coherent states 
it holds AX = l/y/2. 



C. Example: Two-mode squeezed vacuum 

To illustrate these two-mode classical bounds let us 
consider that the measured state is a two-mode squeezed 
vacuum, that in the photon-number basis reads, 



|C) = v / W I Ec"NiN2, 



(5.15) 



n=0 



where we have assumed real parameter £ without loss of 
generality. The statistic of the joint number p n>n and the 
total number p2 n are 



Pi, 



P2r 



(i - c 2 ) c 2 



(5.16) 



while the statistics of the quadrature-difference X = 
X\ — Xi is Gaussian with 



(AX) 



i-C 

2(1 + 0' 



(5.17) 



Numerically we have found that the classical bound 
on joint-number measurements is always violated p n , n > 
p\ n for some n when £ > 0.41. The classical bound for 
total number is never violated since p2n — Pn,n < Pb.in 
for all C- 

More specifically, for n — 1 we havepi.i = (1 — C 2 )C 2 > 
exp(-2) = 0.135 for all C in the interval 0.41 < C < 
0.91. The maximum violation occurs for ( 2 = 1/2 so that 
Pi : i = 0.25, and there is an 85 % violation of the classical 
bound. For the total number we have that the classical 
bound is not surpassed since pb,i = 2exp(— 2) = 0.27. 

On the other hand the classical bound on quadrature 
difference is always surpassed since the statistics is Gaus- 
sian and AX < l/\/2 for all C- For example, for ( 2 = 1/2 
(this is mean total number of photons (n) = 2) we have 
Po = 1.36, while the classical upper bound in Eq. (|5.14p 
is pb = 0.56, so that there is a percentage of violation 
100(po — Pb)/Pb — 141% approximately. The outputs 
x that contravene Eq. (|5. 14|) are all x in the interval 
—0.39 < x < 0.39, which represent an 82% probability 
since J^ 9 39 p x dx = 0.82. 



The above methods can be adapted to situations de- 
scribed by finite-dimensional Hilbert spaces, exemplified 
by spin-j systems. This can be readily done in terms 
of SU(2) Q and P functions, which are defined after the 
SU(2) coherent states \j, CI) as JT3] 

P = J d 2 np(n)\j,n)(j,n\, q(o) = ^±l(j,n| P |i,n), 

(6.1) 

with d 2 Vl — sin 9d6d(f>, and 



U,n>= E 



2 ? . V sin^ m (| l(<- 



j+m ^ 



■exp[-i(j + m)<f>]\j,m), (6.2) 

where |j, m) are the eigenstates of the spin component j'3 
with eigenvalue m, while it > 9 > and tt > <fi > — tt. 
The analog of Eq. (L3]) is 



Pr, 



47T 

2j + l 

47T 
27+T 



d 2 np m (n)Q(n) 
d 2 np{n)Q m (ft), (6.3) 



leading to the following bounds for classical measure- 
ments 



Pr, 



< 



-in 



2.7 + 1 



(6.4) 



(for finite-dimensional systems trA m is always finite), 
while the upper bound for classical states is 



47T „ 
Pm<—Q 



m,max — Pb.7n- 



(6.5) 



By construction the SU(2) coherent states are classical 
both as measured and measuring states. For the POVM 



A n = ^l|i,n)(j,n| 1 trAn = 2 - y _ ' 



47T 



47T 



(6.6) 



the upper bound (|6.4p is satisfied for all measured states 
since 

An 

Pn = Q(O) < 77— — -QmaxtrAo = Q max , (6.7) 
2j + 1 

where Q(tt) is the SU(2) Q function of the measured 
state. Likewise, when the measured state is coherent 
p = \ j, O) (j, fi| we have 



4tt 



Pn 



2j + 1 



Qm(Q), 



(6.8) 



where <5 m (Sl) is the SU(2) Q function of A m , so that 
(|6.5p is satisfied by all measurements. 

For the sake of illustration next we consider some ex- 
amples for the simplest cases j = 1/2, and j = 1. 
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A. j = 1/2 



Nonclassical states 



For j = 1/2 the SU(2) coherent states read, in the 
\j, m) basis, as 



|l/2,fi) = sin(-)|l/2,-l/2) 



- ros ( - ]exp(-#)|l/2,l/2). 



Every state and POVM are of the form 

p=l(ao + r-<T), Q(fi) = i-(l + r-n) : 



(6.10) 

A m = A m (cr + n n ■ tr) , Qm(H) = ^ (1 + r m ■ SI) , 

Zn 

. (6 - H) 

where u are the Pauli matrices in the \j, m) basis, <7o is 
the identity, r, and r m are three-dimensional real vectors 
with |r|, \r m \ < 1, £7 = (sin 9 cos 4>, sin 9 sin 0, cos 6>), and 
A ro > 0. 

Let us show that for j = 1/2 no state exhibits nonclas- 
sical properties. This is because for every A m , and using 



that tr {<Jj(Jk) — 2<5 



tr (pA m ) = A m (1 + r, n 

Z7T 



(6.12) 



which is 



so that p m > pb.m would imply r m ■ r > \r v 
not possible since \r\ < 1. 

Similarly, there are no measurements exceeding the 
classical bounds since p m — A m (l + r m ■ r), Q max = 
(1 + |r|)/(47r), and trA m = 2A m so that the viola- 
tion of the classical bound (|6.4p would be equivalent to 
I'm ■ r > |r|, which is not possible since \r m \ < 1. 

This lack of nonclassical states agrees with the ap- 
proach in Ref. 18] and with the fact that for j = 1/2 
all pure states are SU(2) coherent states. On the other 
hand, this is in sharp contrast with the fact that all pure 
state have negative values of the SU(2) Wigner func- 
tion [l||. Nevertheless, it is worth pointing out that for 
finite-dimensional systems the SU(2) distributions such 
as P{£1) and the Wigner function are nor uniquely de- 
fined in contrast with their infinite-dimensional counter- 
parts [H,[l9|. 



B. j = l 

For j — I the coherent states (|6.2[) arc 
|l,fi) = sin 2 (£)|l,-l) 

+ \/2sinQ^cos^0expH0)|l,O) 



cos^ ( -)expH20)|l,l) 



(6.13) 



As a measurement revealing nonclassical states we can 
consider the projection on the state |1,0) (in the \j, m) 
basis), 



A =p Q = |1,0)<1,0|, 



(6.9) with Q function, 



Q (Q) = — sin 2 ( 



Qo.max n ) 
87T 



(6.14) 



(6.15) 



so that the classical upper bound in Eq. (|6.5p is pbfi = 
1/2. This is clearly violated when the measured state is 
the same state p = po, since in such a case po = 1 > 
Pb,o = 1/2, and there is a 100% violation of the classical 
bound. 

This agrees with the fact that the state |1,0) in Eq. 
(|6.14p can be regarded as the limit of SU(2) squeezed 
states [13, This also agrees with the result in Ref. 
[HI stating that for j = 1 classical behavior (i. e., 
nonsingular positive f (O)) is equivalent to non- negative 
covariance-like matrix 



Zk,e = ((jkji +jdk)) - h,e - (jk)(h)- 

For the state |1, 0) in Eq. we have [2l| 



Z = 




(6.16) 



(6.17) 



so that the state is nonclassical. 



2. Nonclassical measurements 

We can provide an example of nonclassical measure- 
ment. To this end we consider the measurement (in the 
\j,m) basis) A = po = |1,0)(1,0|. As measured state 
we consider the phase averaged equatorial SU(2) coher- 
ent state 

p=^-l #|l,0 = 7r/2,0)<l,0 = 7r/2,</>|, (6.18) 

where \j, 9, <fi) are the SU(2) coherent states. In the 
\j, m)(j, m\ basis this is 

P = \\h 1| + ||1,0><1,0] + -1|, (6.19) 

with Q function 



^^ii( 1 + r in2 ' 



?max - (6.20) 



The probability is po = 1/2, which is 167 % above the 
classical bound in Eq. (|6.4[) , 



-in 



Pb 



2j + l 



ttrAo = -. 



(6.21) 
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VII. CONCLUSIONS 

We have provided feasible practical procedures to re- 
veal the quantum nature of states and measurements. We 
have illustrated them with the most practical measuring 
schemes available, such as photon-number and quadra- 
ture measurements. 

The nonclassical tests proposed in this approach are 
exceedingly simple since definite conclusions are obtained 
without evaluation of moments, or any other more sophis- 
ticated data analysis. This is reflected on the robustness 
of these nonclassical criteria under practical imperfec- 
tions, such as finite detection efficiencies and finite sam- 
pling. 

We have demonstrated that these nonclassical tests are 
independent of other typical quantum signatures such as 
sub-Poissonian statistics, quadrature squeezing, or oscil- 
latory statistics. 
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APPENDIX A: TRACE MEASUREMENTS 

As we have shown above, the classical bounds for 
measurements depend on the trace of the correspond- 
ing POVM elements A m . Incidentally, the trace trA m is 
proportional to the probability of the outcome m when 
the input state is the one of maximum ignorance, p oc 7, 
where 7 is the identity; 



p m (p oc 7) cx trA 

r. 



(Al) 



When trA m depends on m this can be regarded as a kind 
of prior bias, since some outcomes are more probable than 
others even when the measured state presents in principle 
no preference for any outcome. 

In this appendix we present two simple procedures that 
allow us to determine trA m in practice. For simplicity 
and without loss of generality we focus on the single- 
mode case. 

To this end we note that the identity can be expressed 



drrp(r), 



(A2) 



where p(r) are phase-averaged coherent states (which 
correspond to a laser output well above threshold [![) 



p(r) 



1 

2^ 



\r exp(icf))) (r exp(i<p) | , 



(A3) 



and \rexp(i(j>)) are coherent states \a) with a = 
rexp(i(f>). When illuminating the detector with phase- 
averaged coherent states p(r) the measured statistics is 
essentially the phase average of the Q function of A TO , 



Pm(r) = - [ d(f>Q m [r exp(i^)] 



(A4) 



After Eq. (|A2[) the desired trA m can be determined by 
repeating the measurement for different input states p(r) 
by suitably varying the coherent amplitude r 



trA m = tr (A m 7) = 2 / rdrp m (r) 



(A5) 



Although this recall a tomographic reconstruction of A m 
8] , this is not the case since the illuminating state already 
carries the angular integration. 

The same goal can be achieved following an slightly 
different strategy. This is by illuminating the detection 
system with thermal-chaotic states (|2.15j) . The corre- 
sponding Q function (|2.17[) is a Gaussian centered at the 
origin of the complex plane and its width increases when 
the average mean number of photons n tc increases. Thus 
for n tc large enough the Q function of the thermal-chaotic 
state will be approximately constant Q(a) ~ Qq on the 
area where P m (&) ^ 0, so that the statistics will be pro- 
portional to trA m : 

Pm = Tr J d 2 aP m (a)Q(a) 

~ nQo J d 2 aP m (a) = 7rQ trA m , (A6) 



and 



trA, 



Pn 



lim 

n ta — »oo TTLJq 



(A7) 



To illustrate this idea let us consider a one-photon de- 
tector whose output is contaminated by the vacuum and 
two-photon contributions 



Ai = q 



+ p|l)(l|+flf|2)<2|, 



(A8) 



with trAi = p + 2q. When illuminated with the thermal- 
chaotic state we get 



1 



Pi 



ntc 



q + p 



n tc 



n tc + 1 (n tc + I) 2 



(A9) 



For example, for p = 1, and q = 0.1, and taking ttQq 
1/ (n tc + 1), we get for n tc = 100 



trAi = 1.2, 
i. e., only 1 % error. 



Pi 



= 1.19, 



(A10) 
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